Dynamic molecular mechanism of the nuclear pore complex permeability barrier.

Nuclear pore complexes (NPCs) mediate nucleocytoplasmic transport of specific macromolecules while impeding the exchange of unsolicited material. However, key aspects of this gating mechanism remain controversial. To address this issue, we determined the nanoscopic behavior of the permeability barrier directly within yeast S. cerevisiae NPCs at transport-relevant timescales. We show that the large intrinsically disordered domains of phenylalanine-glycine repeat nucleoporins (FG Nups) exhibit highly dynamic fluctuations to create transient voids in the permeability barrier that continuously shape-shift and reseal, resembling a radial polymer brush. Together with cargo-carrying transport factors the FG domains form a feature called the central plug, which is also highly dynamic. Remarkably, NPC mutants with longer FG domains show interweaving meshwork-like behavior that attenuates nucleocytoplasmic transport in vivo. Importantly, the bona fide nanoscale NPC behaviors and morphologies are not recapitulated by in vitro FG domain hydrogels. NPCs also exclude self-assembling FG domain condensates in vivo, thereby indicating that the permeability barrier is not generated by a self-assembling phase condensate, but rather is largely a polymer brush, organized by the NPC scaffold, whose dynamic gating selectivity is strongly enhanced by the presence of transport factors.

Ultrafast Transverse Modulation of Free Electrons by Interaction with Shaped Optical Fields.

Spatiotemporal electron-beam shaping is a bold frontier of electron microscopy. Over the past decade, shaping methods evolved from static phase plates to low-speed electrostatic and magnetostatic displays. Recently, a swift change of paradigm utilizing light to control free electrons has emerged. Here, we experimentally demonstrate arbitrary transverse modulation of electron beams without complicated electron-optics elements or material nanostructures, but rather using shaped light beams. On-demand spatial modulation of electron wavepackets is obtained via inelastic interaction with transversely shaped ultrafast light fields controlled by an external spatial light modulator. We illustrate this method for the cases of Hermite-Gaussian and Laguerre-Gaussian modulation and discuss their use in enhancing microscope sensitivity. Our approach dramatically widens the range of patterns that can be imprinted on the electron profile and greatly facilitates tailored electron-beam shaping.

Motional clustering in supra-τc conformational exchange influences NOE cross-relaxation rate.

Biomolecular spin relaxation processes, such as the NOE, are commonly modeled by rotational τc-tumbling combined with fast motions on the sub-τc timescale. Motions on the supra-τc timescale, in contrast, are considered to be completely decorrelated to the molecular tumbling and therefore invisible. Here, we show how supra-τc dynamics can nonetheless influence the NOE build-up between methyl groups. This effect arises because supra-τc motions can cluster the fast-motion ensembles into discrete states, affecting distance averaging as well as the fast-motion order parameter and hence the cross-relaxation rate. We present a computational approach to estimate methyl-methyl cross-relaxation rates from extensive (>100×τc) all-atom molecular dynamics (MD) trajectories on the example of the 723-residue protein Malate Synthase G. The approach uses Markov state models (MSMs) to resolve transitions between metastable states and thus to discriminate between sub-τc and supra-τc conformational exchange. We find that supra-τc exchange typically increases NOESY cross-peak intensities. The methods described in this work extend the theory of modeling sub-µs dynamics in spin relaxation and thus contribute to a quantitative estimation of NOE cross-relaxation rates from MD simulations, eventually leading to increased precision in structural and functional studies of large proteins.

Insights into SusCD-mediated glycan import by a prominent gut symbiont.

In Bacteroidetes, one of the dominant phyla of the mammalian gut, active uptake of large nutrients across the outer membrane is mediated by SusCD protein complexes via a "pedal bin" transport mechanism. However, many features of SusCD function in glycan uptake remain unclear, including ligand binding, the role of the SusD lid and the size limit for substrate transport. Here we characterise the ß2,6 fructo-oligosaccharide (FOS) importing SusCD from Bacteroides thetaiotaomicron (Bt1762-Bt1763) to shed light on SusCD function. Co-crystal structures reveal residues involved in glycan recognition and suggest that the large binding cavity can accommodate several substrate molecules, each up to ~2.5 kDa in size, a finding supported by native mass spectrometry and isothermal titration calorimetry. Mutational studies in vivo provide functional insights into the key structural features of the SusCD apparatus and cryo-EM of the intact dimeric SusCD complex reveals several distinct states of the transporter, directly visualising the dynamics of the pedal bin transport mechanism.

Enhancing NMR derived ensembles with kinetics on multiple timescales.

Nuclear magnetic resonance (NMR) has the unique advantage of elucidating the structure and dynamics of biomolecules in solution at physiological temperatures, where they are in constant movement on timescales from picoseconds to milliseconds. Such motions have been shown to be critical for enzyme catalysis, allosteric regulation, and molecular recognition. With NMR being particularly sensitive to these timescales, detailed information about the kinetics can be acquired. However, nearly all methods of NMR-based biomolecular structure determination neglect kinetics, which introduces a large approximation to the underlying physics, limiting both structural resolution and the ability to accurately determine molecular flexibility. Here we present the Kinetic Ensemble approach that uses a hierarchy of interconversion rates between a set of ensemble members to rigorously calculate Nuclear Overhauser Effect (NOE) intensities. It can be used to simultaneously refine both temporal and structural coordinates. By generalizing ideas from the extended model free approach, the method can analyze the amplitudes and kinetics of motions anywhere along the backbone or side chains. Furthermore, analysis of a large set of crystal structures suggests that NOE data contains a surprising amount of high-resolution information that is better modeled using our approach. The Kinetic Ensemble approach provides the means to unify numerous types of experiments under a single quantitative framework and more fully characterize and exploit kinetically distinct protein states. While we apply the approach here to the protein ubiquitin and cross validate it with previously derived datasets, the approach can be applied to any protein for which NOE data is available.

Regulation of α-synuclein by chaperones in mammalian cells.

Neurodegeneration in patients with Parkinson's disease is correlated with the occurrence of Lewy bodies-intracellular inclusions that contain aggregates of the intrinsically disordered protein α-synuclein1. The aggregation propensity of α-synuclein in cells is modulated by specific factors that include post-translational modifications2,3, Abelson-kinase-mediated phosphorylation4,5 and interactions with intracellular machineries such as molecular chaperones, although the underlying mechanisms are unclear6-8. Here we systematically characterize the interaction of molecular chaperones with α-synuclein in vitro as well as in cells at the atomic level. We find that six highly divergent molecular chaperones commonly recognize a canonical motif in α-synuclein, consisting of the N terminus and a segment around Tyr39, and hinder the aggregation of α-synuclein. NMR experiments9 in cells show that the same transient interaction pattern is preserved inside living mammalian cells. Specific inhibition of the interactions between α-synuclein and the chaperone HSC70 and members of the HSP90 family, including HSP90ß, results in transient membrane binding and triggers a remarkable re-localization of α-synuclein to the mitochondria and concomitant formation of aggregates. Phosphorylation of α-synuclein at Tyr39 directly impairs the interaction of α-synuclein with chaperones, thus providing a functional explanation for the role of Abelson kinase in Parkinson's disease. Our results establish a master regulatory mechanism of α-synuclein function and aggregation in mammalian cells, extending the functional repertoire of molecular chaperones and highlighting new perspectives for therapeutic interventions for Parkinson's disease.

An integrative protocol for the structure determination of the mouse ASC-PYD filament.

Several inflammasomes that activate as part of the eukaryotic innate immune response contain long helical filaments of the adaptor protein ASC as a central structural element. Here, we describe a technical protocol that has enabled the structure determination of the filament of the ASC pyrin domain. The protocol integrates data from cryo-electron microscopy and solid-state NMR spectroscopy into a single simulated annealing protocol to determine structural coordinates that fit all input data optimally. The structure shows that the ASC pyrin domain filament is formed by helical stacking of individual pyrin domains forms and that the CARD domains are flexibly attached to the filament outside. An artificial perturbation of the input data shows that the integrated structure determination protocol can allow high quality structures even at resolutions of the electron density map as low 8Å. The protocol is extendable to other structural input data from biochemical or biophysical experiments.

The dynamic dimer structure of the chaperone Trigger Factor.

The chaperone Trigger Factor (TF) from Escherichia coli forms a dimer at cellular concentrations. While the monomer structure of TF is well known, the spatial arrangement of this dimeric chaperone storage form has remained unclear. Here, we determine its structure by a combination of high-resolution NMR spectroscopy and biophysical methods. TF forms a symmetric head-to-tail dimer, where the ribosome binding domain is in contact with the substrate binding domain, while the peptidyl-prolyl isomerase domain contributes only slightly to the dimer affinity. The dimer structure is highly dynamic, with the two ribosome binding domains populating a conformational ensemble in the center. These dynamics result from intermolecular in trans interactions of the TF client-binding site with the ribosome binding domain, which is conformationally frustrated in the absence of the ribosome. The avidity in the dimer structure explains how the dimeric state of TF can be monomerized also by weakly interacting clients.

High-Speed Atomic Force Microscopy Visualization of the Dynamics of the Multienzyme Fatty Acid Synthase.

Multienzymes, such as the protein metazoan fatty acid synthase (FAS), are giant and highly dynamic molecular machines for critical biosynthetic processes. The molecular architecture of FAS was elucidated by static high-resolution crystallographic analysis, while electron microscopy revealed large-scale conformational variability in FAS with some correlation to functional states in catalysis. However, little is known about time scales of conformational dynamics, the trajectory of motions in individual FAS molecules, and the extent of coupling between catalysis and structural changes. Here, we present an experimental single-molecule approach to film immobilized or selectively tethered FAS in solution at different viewing angles and high spatiotemporal resolution using high-speed atomic force microscopy. Mobility of individual regions of the multienzyme is recognized in video sequences, and correlation of shape features implies a convergence of temporal resolution and velocity of FAS dynamics. Conformational variety can be identified and grouped by reference-free 2D class averaging, enabling the tracking of conformational transitions in movies. The approach presented here is suited for comprehensive studies of the dynamics of FAS and other multienzymes in aqueous solution at the single-molecule level.

The Pharmacological Profile of a Novel Highly Potent Bisphosphonate, OX14 (1-Fluoro-2-(Imidazo-[1,2-α]Pyridin-3-yl)-Ethyl-Bisphosphonate).

Bisphosphonates are widely used in the treatment of clinical disorders characterized by increased bone resorption, including osteoporosis, Paget's disease, and the skeletal complications of malignancy. The antiresorptive potency of the nitrogen-containing bisphosphonates on bone in vivo is now recognized to depend upon two key properties, namely mineral binding affinity and inhibitory activity on farnesyl pyrophosphate synthase (FPPS), and these properties vary independently of each other in individual bisphosphonates. The better understanding of structure activity relationships among the bisphosphonates has enabled us to design a series of novel bisphosphonates with a range of mineral binding properties and antiresorptive potencies. Among these is a highly potent bisphosphonate, 1-fluoro-2-(imidazo-[1,2 alpha]pyridin-3-yl)-ethyl-bisphosphonate, also known as OX14, which is a strong inhibitor of FPPS, but has lower binding affinity for bone mineral than most of the commonly studied bisphosphonates. The aim of this work was to characterize OX14 pharmacologically in relation to several of the bisphosphonates currently used clinically. When OX14 was compared to zoledronate (ZOL), risedronate (RIS), and minodronate (MIN), it was as potent at inhibiting FPPS in vitro but had significantly lower binding affinity to hydroxyapatite (HAP) columns than ALN, ZOL, RIS, and MIN. When injected i.v. into growing Sprague Dawley rats, OX14 was excreted into the urine to a greater extent than the other bisphosphonates, indicating reduced short-term skeletal uptake and retention. In studies in both Sprague Dawley rats and C57BL/6J mice, OX14 inhibited bone resorption, with an antiresorptive potency equivalent to or greater than the comparator bisphosphonates. In the JJN3-NSG murine model of myeloma-induced bone disease, OX14 significantly prevented the formation of osteolytic lesions (p < 0.05). In summary, OX14 is a new, highly potent bisphosphonate with lower bone binding affinity than other clinically relevant bisphosphonates. This renders OX14 an interesting potential candidate for further development for its potential skeletal and nonskeletal benefits. © 2017 American Society for Bone and Mineral Research.

Spatiotemporal dynamics of the nuclear pore complex transport barrier resolved by high-speed atomic force microscopy.

Nuclear pore complexes (NPCs) are biological nanomachines that mediate the bidirectional traffic of macromolecules between the cytoplasm and nucleus in eukaryotic cells. This process involves numerous intrinsically disordered, barrier-forming proteins known as phenylalanine-glycine nucleoporins (FG Nups) that are tethered inside each pore. The selective barrier mechanism has so far remained unresolved because the FG Nups have eluded direct structural analysis within NPCs. Here, high-speed atomic force microscopy is used to visualize the nanoscopic spatiotemporal dynamics of FG Nups inside Xenopus laevis oocyte NPCs at timescales of ∼100 ms. Our results show that the cytoplasmic orifice is circumscribed by highly flexible, dynamically fluctuating FG Nups that rapidly elongate and retract, consistent with the diffusive motion of tethered polypeptide chains. On this basis, intermingling FG Nups exhibit transient entanglements in the central channel, but do not cohere into a tightly crosslinked meshwork. Therefore, the basic functional form of the NPC barrier is comprised of highly dynamic FG Nups that manifest as a central plug or transporter when averaged in space and time.

Intrinsic regulation of FIC-domain AMP-transferases by oligomerization and automodification.

Filamentation induced by cyclic AMP (FIC)-domain enzymes catalyze adenylylation or other posttranslational modifications of target proteins to control their function. Recently, we have shown that Fic enzymes are autoinhibited by an α-helix (αinh) that partly obstructs the active site. For the single-domain class III Fic proteins, the αinh is located at the C terminus and its deletion relieves autoinhibition. However, it has remained unclear how activation occurs naturally. Here, we show by structural, biophysical, and enzymatic analyses combined with in vivo data that the class III Fic protein NmFic from Neisseria meningitidis gets autoadenylylated in cis, thereby autonomously relieving autoinhibition and thus allowing subsequent adenylylation of its target, the DNA gyrase subunit GyrB. Furthermore, we show that NmFic activation is antagonized by tetramerization. The combination of autoadenylylation and tetramerization results in nonmonotonic concentration dependence of NmFic activity and a pronounced lag phase in the progress of target adenylylation. Bioinformatic analyses indicate that this elaborate dual-control mechanism is conserved throughout class III Fic proteins.

A molecular mechanism of chaperone-client recognition.

Molecular chaperones are essential in aiding client proteins to fold into their native structure and in maintaining cellular protein homeostasis. However, mechanistic aspects of chaperone function are still not well understood at the atomic level. We use nuclear magnetic resonance spectroscopy to elucidate the mechanism underlying client recognition by the adenosine triphosphate-independent chaperone Spy at the atomic level and derive a structural model for the chaperone-client complex. Spy interacts with its partially folded client Im7 by selective recognition of flexible, locally frustrated regions in a dynamic fashion. The interaction with Spy destabilizes a partially folded client but spatially compacts an unfolded client conformational ensemble. By increasing client backbone dynamics, the chaperone facilitates the search for the native structure. A comparison of the interaction of Im7 with two other chaperones suggests that the underlying principle of recognizing frustrated segments is of a fundamental nature.

Sequence-specific solid-state NMR assignments of the mouse ASC PYRIN domain in its filament form.

The apoptosis-associated speck-like protein (ASC protein) plays a central role in eukaryotic innate immune response. Upon infection, multiple ASC molecules assemble into long filaments, which are fundamental for triggering the cellular defense mechanism by starting an inflammatory cascade with the activation of caspase-1. ASC is composed of two domains, the C-terminal caspase-recruitment domain, which is involved in the recruitment of the caspase, and the N-terminal PYRIN domain (PYD), which is responsible for the formation of the filament. Here we present the (13)C and (15)N chemical shift assignment for filaments formed by the PYD of mouse ASC, a 91-residue protein. The backbone between residues 4 and 84 is assigned without interruption. Also, 86 % of the sidechain resonances for this stretch are assigned. Residues 1-3 and 85-91 show unfavorable dynamics and are not observed. Secondary chemical-shift analysis shows the presence of six α-helices.

Expression and Genetic Activation of Cyclic Di-GMP-Specific Phosphodiesterases in Escherichia coli.

UNLABELLED: Intracellular levels of the bacterial second messenger cyclic di-GMP (c-di-GMP) are controlled by antagonistic activities of diguanylate cyclases and phosphodiesterases. The phosphodiesterase PdeH was identified as a key regulator of motility in Escherichia coli, while deletions of any of the other 12 genes encoding potential phosphodiesterases did not interfere with motility. To analyze the roles of E. coli phosphodiesterases, we demonstrated that most of these proteins are expressed under laboratory conditions. We next isolated suppressor mutations in six phosphodiesterase genes, which reinstate motility in the absence of PdeH by reducing cellular levels of c-di-GMP. Expression of all mutant alleles also led to a reduction of biofilm formation. Thus, all of these proteins are bona fide phosphodiesterases that are capable of interfering with different c-di-GMP-responsive output systems by affecting the global c-di-GMP pool. This argues that E. coli possesses several phosphodiesterases that are inactive under laboratory conditions because they lack appropriate input signals. Finally, one of these phosphodiesterases, PdeL, was studied in more detail. We demonstrated that this protein acts as a transcription factor to control its own expression. Motile suppressor alleles led to a strong increase of PdeL activity and elevated pdeL transcription, suggesting that enzymatic activity and transcriptional control are coupled. In agreement with this, we showed that overall cellular levels of c-di-GMP control pdeL transcription and that this control depends on PdeL itself. We thus propose that PdeL acts both as an enzyme and as a c-di-GMP sensor to couple transcriptional activity to the c-di-GMP status of the cell. IMPORTANCE: Most bacteria possess multiple diguanylate cyclases and phosphodiesterases. Genetic studies have proposed that these enzymes show signaling specificity by contributing to distinct cellular processes without much cross talk. Thus, spatial separation of individual c-di-GMP signaling units was postulated. However, since most cyclases and phosphodiesterases harbor N-terminal signal input domains, it is equally possible that most of these enzymes lack their activating signals under laboratory conditions, thereby simulating signaling specificity on a genetic level. We demonstrate that a subset of E. coli phosphodiesterases can be activated genetically to affect the global c-di-GMP pool and thus influence different c-di-GMP-dependent processes. Although this does not exclude spatial confinement of individual phosphodiesterases, this study emphasizes the importance of environmental signals for activation of phosphodiesterases.

A Rigorous and Efficient Method To Reweight Very Large Conformational Ensembles Using Average Experimental Data and To Determine Their Relative Information Content.

Flexible polypeptides such as unfolded proteins may access an astronomical number of conformations. The most advanced simulations of such states usually comprise tens of thousands of individual structures. In principle, a comparison of parameters predicted from such ensembles to experimental data provides a measure of their quality. In practice, analyses that go beyond the comparison of unbiased average data have been impossible to carry out on the entirety of such very large ensembles and have, therefore, been restricted to much smaller subensembles and/or nondeterministic algorithms. Here, we show that such very large ensembles, on the order of 10(4) to 10(5) conformations, can be analyzed in full by a maximum entropy fit to experimental average data. Maximizing the entropy of the population weights of individual conformations under experimental χ(2) constraints is a convex optimization problem, which can be solved in a very efficient and robust manner to a unique global solution even for very large ensembles. Since the population weights can be determined reliably, the reweighted full ensemble presents the best model of the combined information from simulation and experiment. Furthermore, since the reduction of entropy due to the experimental constraints is well-defined, its value provides a robust measure of the information content of the experimental data relative to the simulated ensemble and an indication for the density of the sampling of conformational space. The method is applied to the reweighting of a 35,000 frame molecular dynamics trajectory of the nonapeptide EGAAWAASS by extensive NMR (3)J coupling and RDC data. The analysis shows that RDCs provide significantly more information than (3)J couplings and that a discontinuity in the RDC pattern at the central tryptophan is caused by a cluster of helical conformations. Reweighting factors are moderate and consistent with errors in MD force fields of less than 3kT. The required reweighting is larger for an ensemble derived from a statistical coil model, consistent with its coarser nature. We call the method COPER, for convex optimization for ensemble reweighting. Similar advantages of large-scale efficiency and robustness can be obtained for other ensemble analysis methods with convex targets and constraints, such as constrained χ(2) minimization and the maximum occurrence method.

Structure and assembly of the mouse ASC inflammasome by combined NMR spectroscopy and cryo-electron microscopy.

Inflammasomes are multiprotein complexes that control the innate immune response by activating caspase-1, thus promoting the secretion of cytokines in response to invading pathogens and endogenous triggers. Assembly of inflammasomes is induced by activation of a receptor protein. Many inflammasome receptors require the adapter protein ASC [apoptosis-associated speck-like protein containing a caspase-recruitment domain (CARD)], which consists of two domains, the N-terminal pyrin domain (PYD) and the C-terminal CARD. Upon activation, ASC forms large oligomeric filaments, which facilitate procaspase-1 recruitment. Here, we characterize the structure and filament formation of mouse ASC in vitro at atomic resolution. Information from cryo-electron microscopy and solid-state NMR spectroscopy is combined in a single structure calculation to obtain the atomic-resolution structure of the ASC filament. Perturbations of NMR resonances upon filament formation monitor the specific binding interfaces of ASC-PYD association. Importantly, NMR experiments show the rigidity of the PYD forming the core of the filament as well as the high mobility of the CARD relative to this core. The findings are validated by structure-based mutagenesis experiments in cultured macrophages. The 3D structure of the mouse ASC-PYD filament is highly similar to the recently determined human ASC-PYD filament, suggesting evolutionary conservation of ASC-dependent inflammasome mechanisms.

On-the-Fly Integration of Data from a Spin-Diffusion-Based NMR Experiment into Protein-Ligand Docking.

INPHARMA (interligand nuclear Overhauser enhancement for pharmacophore mapping) determines the relative orientation of two competitive ligands in the protein binding pocket. It is based on the observation of interligand transferred NOEs mediated by spin diffusion through protons of the protein and is, therefore, sensitive to the specific interactions of each of the two ligands with the protein. We show how this information can be directly included into a protein-ligand docking program to guide the prediction of the complex structures. Agreement between the experimental and back-calculated spectra based on the full relaxation matrix approach is translated into a score contribution that is combined with the scoring function ChemPLP of our docking tool PLANTS. This combined score is then used to predict the poses of five weakly bound cAMP-dependent protein kinase (PKA) ligands. After optimizing the setup, which finally also included trNOE data and optimized protonation states, very good success rates were obtained for all combinations of three ligands. For one additional ligand, no conclusive results could be obtained due to the ambiguous electron density of the ligand in the X-ray structure, which does not disprove alternative ligand poses. The failures of the remaining ligand are caused by suboptimal locations of specific protein side chains. Therefore, side-chain flexibility should be included in an improved INPHARMA-PLANTS version. This will reduce the strong dependence on the used protein input structure leading to improved scores overall, not only for this last ligand.

A combination of spin diffusion methods for the determination of protein-ligand complex structural ensembles.

Structure-based drug design (SBDD) is a powerful and widely used approach to optimize affinity of drug candidates. With the recently introduced INPHARMA method, the binding mode of small molecules to their protein target can be characterized even if no spectroscopic information about the protein is known. Here, we show that the combination of the spin-diffusion-based NMR methods INPHARMA, trNOE, and STD results in an accurate scoring function for docking modes and therefore determination of protein-ligand complex structures. Applications are shown on the model system protein kinase A and the drug targets glycogen phosphorylase and soluble epoxide hydrolase (sEH). Multiplexing of several ligands improves the reliability of the scoring function further. The new score allows in the case of sEH detecting two binding modes of the ligand in its binding site, which was corroborated by X-ray analysis.

Sampling of Glycan-Bound Conformers by the Anti-HIV Lectin Oscillatoria agardhii agglutinin in the Absence of Sugar.

Lectins from different sources have been shown to interfere with HIV infection by binding to the sugars of viral-envelope glycoproteins. Three-dimensional atomic structures of a number of HIV-inactivating lectins have been determined, both as free proteins and in glycan-bound forms. However, details on the mechanism of recognition and binding to sugars are elusive. Herein we focus on the anti-HIV lectin OAA from Oscillatoria agardhii: We show that in the absence of sugars in solution, both the sugar-free and sugar-bound protein conformations that were observed in the X-ray crystal structures exist as conformational substates. Our results suggest that glycan recognition occurs by conformational selection within the ground state; this model differs from the popular "excited-state" model. Our findings provide further insight into molecular recognition of the major receptor on the HIV virus by OAA. These details can potentially be used for the optimization and/or development of preventive anti-HIV therapeutics.